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0.-  FOREWORD 


The  work  done  during  the  period  covered  by  the  present  FINAL  Scientific 
Report  pertains  to  several  topics  related  to  the  application  of  functio- 
nal analysis  in  fluid-mechanics. 

It  has  lead  to  the  following  papers  : 

1. -  L.G.  NAPOLITANO  : "Functional  Analysis  Approach  for  the  Derivation 

of  Hybrid  Variational  Functionals". 

This  paper  was  presented  at  the  IUTAM/IUM  Symposium  on  Applications 
of  Methods  of  Functional  Analysis  to  Problems  of  Mechanics  held  in 
Marseille,  France,  on  September  1975. 

The  Proceedings  have  already  been  published  by  Springer  Verlag.  The 
complete  reference  is  : 

"Applications  of  Methods  of  Functional  Analysis  to  Problems  in  Mechanics" 
Edited  by  P.  Germain  and  B.  Nayroles;  Lecture  Notes  in  Mathematics,  no. 
503  , Springer-Verlag,  Berlin  - Heidelberg  - New  York,  1976. 

In  the  paper  acknowledgment  to  AFOSR  sponsorship  was  incorrectly  stated 
as  pertaining  to  Grant  No.  AFOSR  74-2704.  The  correct  statement  should 
have  been  Grant  No.  AF0SR-76-2889. 

2. -  L.G.  NAPOLITANO  :"The  functional  inverse  approach  for  the  variational 

formulation  of  boundary  value  problems". 

This  paper  was  presented  at  the  Congress  on  Modern  Problems  in  Continuum 
Mechanics  held  in  Turin  (Italy)  on  November  1975. 


The  proceedings  of  this  Meeting  are  in  press. 


3. -  L.G.  NAPOLITANO  :"Functional  Analysis  Derivation  of  Hybrid  Variationa 

Functionals  for  Fourth  Order  Elliptical  Operators". 

This  paper  was  presented  at  the  Second  International  Symposium  on  Finite 
Elements  Methods  in  Flow  Problems  held  in  S.  Margherita  Ligure  (Italy) 
on  June  1976. 

The  Preprints  of  the  Symposium  have  been  published  already  and  copy  of 
them  can  be  obtained  from  ICCAD  (International  Centre  for  Computer  Ai- 
ded Design)  Genova  (Italy). 

4. -  L.G.  NAPOLITANO  and  V.  LOSITO  :"The  Closed  Spline  Functions". 

This  paper  will  be  submitted  for  publication  to  an  International  Journal 

Since  papers  No.  1 and  3 above  are  readily  available  (and  an  appropriate 
number  of  copies  have  already  been  forwarded  to  the  Sponsoring  Agency) 
the  subject  Final  Report  will  contain,  in  extenso,  only  the  papers  No. 

2 and  4 above  which  constitute  Parts  I and  II  respectively  of  the  Report 
itself. 


PART  I 

THE  FUNCTIONAL  INVERSE  APPROACH  FOR  THE  VARIATIONAL 
FORMULATION  OF  BOUNDARY  VALUE  PROBLEMS 


V 

by 

L.  G.  Napolitano 
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THE  FUNCTIONAL  INVERSE  APPROACH  FOR  THE  VARIATIONAL 


FORMULATION  OF  BOUNDARY  VALUE  PROBLEMS  (°) 
by 

L.G.  Napolitano 


1.  INTRODUCTION 

The  present  author  has  developed  a functional  approach  for  the  unified 
derivation  of  hybrid  and/or  classical  variational  formulations  of  boun- 
dary value  problems  1 1 ] , C?J.[3Tl<]  ■ 

This  approach, referred  to  as  the  " direct  approach",  considers  as  aivc-n 
a formally  self-adjoint  linear  boundary  value  problem  and  derives  its  dif- 
ferent  variational  formulation  from  the  stationary  properties  of  a twu- 
Field  functional  defined  over  two  suitably  defined  linear  varieties.  The 
direct  approach  hinges  on  an  assumed  factorization  of  the  differential  ope- 
rator defining  the  boundary  value  problem  and,  ultimately,  on  the  corre- 
sponding Gauss-formula. 

A more  powerful  and  general  approach,  to  be  referred  to  as  "inverse  approach 
will  be  presented  here. 

In  dealing  with  a general  variational  theory  of  boundary  value  problems  Lion 
and  Magenes  £5]  have  already  pointed  out  the  connection  existing,  via  an 
underlying  Gauss  formula,  between  the  classes  of  problems  that  can  be  analy- 
sed and  the  assumed  definition  of  an  inner  product. 

The  contention  here  is  that^ ultimately  and  more  generally,  the  connection 


(°)  Paper  presented  at  the  Congress  on  Modern  Problems  in  Continuum 
Mechanics,  Torino,  1976 


• 1 


is  to  be  established  between  an  abstract  Gauss  formula  and  classes  of  boundary 
value  problems.  The  inverse  approach  both  accomplishes  these  tasks  and 
duely  exploits  its  results. 

The  starting  point  is  an  abstract  Gauss  formula  and  the  end  results  are  varia- 
tional formulations  (derived  by  the  method  developed  in  the  direct  approach; 
of  the  classes  of  boundary  value  problems  that  can  be  "associated"  to  the  given 
Gauss  formula. 

As  it  will  be  shown,  an  abstract  Gauss  formula  involves  essentially  two  Hilbert 
spaces  of  elements  If  and  2def ined  over  a domain  a differential  operator 
such  that  if] is  a subset  a number  of  boundary  spaces  defined  over  parti- 

tions of  the  boundary  ofJ2.and  a set  of  principal  and  boundary  operators. 

The  "degrees  of  freedom"  it  affords  thus  pertain  not  only  to  the  choice  of  the 
definition  of  inner  products  but  also  to  the  very  same  choices  of  the  elements  i’f 
cf  the  operator  G and  of  their  tensorial  orders. 

The  power  and  generality  of  the  inverse  approach  stems  from  the  unified  deriva- 
tion of  variational  formulations  for  the  wide  classes  of  boundary  value  problems 
resulting  from  such  diversified  choices. 

The  inverse  approach  presented  in  this  paper  is  a natural  outcome  of  the  di- 
rect approach  previously  developed  by  the  author.  It  thus  makes  uses  of  many 
results  already  established  in  previous  works  to  which  reference  is  made  for 
greater  details  and/or  for  pertinent  discussion  of  the  significance  of  the 
main  assumption  which  underlies  its  development. 


CV 


2.-  FUNDAMENTALS  OF  THE  INVERSE  APPROACH 


2.1  - The  abstract  Gauss  formula 


Let  be  an  open  bounded  subset  of&with  sufficiently  smooth  boundary  O 
and  let  ttxre.  p * Q j ^ nj  C ^ ~ J ~ be  (m  > 1 ) arbitrary 
partitions,  not  necessarily  all  distinct^of  OJ2. 

Let  be  two  appropriate  subsets  of  the  real 

Hilbert  space  L (&J where  the  real  function  if ]/  can  be  of  any 
tensorial  order,  and  G be  a differential  operator  of  order  m and  arbitrary 
tensorial  order.  Assume  that  there  is  a Gauss's  formula  of  the  tvpe  : 

< 6 ^ r > --  ( ■*,<?*)- jL  e;  ( V- ir  A !j  t)  G-  0 

where  (j  V ^ ^ C.  it,  (&J&  subset  of  the  / (^^pace  appropriate  to 
their  common  tensorial  ordery  < , > and  C j _)  denote  inner  pro- 
ducts in  ■Ui  and  \/  , respectively,  G*”is  the  formal  adjoint  of  G and 
the  4 are  bilinear  symmetrical  forms  onTf.^  which  characterize  the  set 
of  (m)  "principal"  (resp.  "natural")  boundary  differential  operators 

P ‘ (resp.  N)  of  order  (j  - 1)  ^resD.  m - j 7 . 

1 J 

The  conditions  on^^  and  its  partitions  necessary  for  the  (formal)  vali- 
dity of  the  Gauss ' formula  are  assumed  to  holdlbT 


The  nature  of  the  elements  if  and  2T  (and  thus  of  the  operator  G)  and 
the  inner  products  -c  / > and  C , )|[and  thus  of  the  sets  of  boundary  diffe- 
rential operators J are  left  "unspecified". 

Denoti«§  by  ^.and  t ri.  the  restrictions  of  £*to\^£nd  *7^,  ■ j[,  respecti 
vely,  and  let  \ ^ ^ 


¥ 


- 2 - 


, \ y a l ?‘V  n:)  = ^tc  ( 

f (?■&  N-z);£-  V '4 

t«- -V \ > i i-‘  i 


where  : 


p -if.,  h, 


n - L*i,  '* 


'Ll 

, *»,  J 


Then  the  Gauss  formula  ( l.  / ) can  be  rewritten  as  : 

(v,6*zj  - if.  " z)  - *-)  tz-Z) 

The  thus  fas  arbitrary  element  2"  6 U,  ( ) is  now  restricted 

to  the  subset  defined  by 

Uz  - I 6-Z.j  ^ j) 

> , { <re  U,  (Ji)  I 6 ° - Vr/sjf  * f 

&■  if  and  (T  will  be  referred^vfespecti vely^  bft  as  the  compatible  (or 
constitutive)  and  non-compatible  parts  of  2?  • 

The  subset  C U,  C%)  supposed  non-empty  [ See [/i] for  relevant  examples  / 
is  orthogonal  to  the  subset  {&vj,  orthogonality  being  understood  with 
respect  to  the  inner  product  <T  y > . 

On  account  of  this,  for  2",  Z 6 the  Gauss  formula  yields  j 

wit k G*6  V~-  * <=  K'(ol)  y h/ut,  6>Lvusl£  KiJL)  : 


L 


In  principle,  the  inner  product  CT  '>  could  be  further  decomposed, 
by  relating  it  to  the  inner  product  pertaining  to  still  another  space 
and  to  other  bilinear  forms.  This  development  is  not  to  be  pursued  here. 

It  will  only  be  supposed  that  there  is  an  element  , belonging  to  the 
cartesian  product  space  , appropriate 

to  the  tensorial  order  of  which  determines  rT  (r  <2.  via  a pertinent 
boundary  value  problem. 


The  elements  *>,  p,  ^ and  will  be  interpreted  as  "generalized"  con- 
straints imposed  on  the  element  2T  of  and  it  will  be  assumed 

that  a given  set  of  values  for  the  generalized  constraints  characterize^ 
a unique  element  (j-) 2. ) as  formalized  in  the  next 

paragraphs . 


2 2/C°nstra''nt  and  solutions  spaces. 


Define 


L h 


n , t 7 


as  a "constraint  element".  Its  components  will  be  referred  to,  respective- 
ly, as  source,  principal  boundary,  natural  boundary  and  constitutive 
(or  compatibility)  constraints.  Thus  they  prescribe,  in  the  order, 
the  source  (s)  in  SL  of  an  element  2r  , its  principal  boundary 

values  p-  (1  - j -m)  on  the  (m)  non  necessarily  coincident  subsets  r0iP . J2, 


of  i) 


; its  (m)  natural  boundary  values  *,‘on  the  (m)  complementary 


subsets  c->t*  Jl  , and,  finally,  the  value  (r)  characterizing  its"non- 

V — - 

compatible"  part  n . 


- 4 - 


It  is  assumed  that  there  exist  suitable  spaces  and  , 

closed  with  respect  to  appropriate  norms,  such  that  if 

f 6D  B jo  fc  Bn.®'  // fc 

then  for  any  C 6 (i.e.  any  ^ t p 6 >7  6 f.) 

the  linear  mapping  F;  ^ ->  llz  C^-J  is  one  to  one;  continuous  and 
bounded . . 

The  space  t>  will  be  called  the  constraint  space  and  its  (closed) 

F-image  H-  - F(&JC  Uz  CJL  J will  be  cal  ed  the  "solution 
space".  F is  an  isomorphism  between  the  constraint  space  ^ and  the  solution 
space  //  . 

The  hypothesis  on  F is  related  to  the  existence  and  uniqueness  of  the  so- 
lution of  the  following  four  linear  boundary  value  problems  : 


< 


> *c  ~ O . 


- 6 - 


when  accounting  for  the  fact  that  r - O 

The  inner  product  among  elements  of  each  one  of  the  four  basic  subspaces 
reduces  to  only  one  of  the  four  contributions  contained  in  the  Gauss 
formula.  Thus,  alternatively  and  equivalently,  the  basic  subspaces  could 
have  been  introduced  directly  by  imposing  such  a requirement  of  the 
Gauss 'formula . 


7.3  The  associated  boundary  value  problem  for  compatible  elements. 

The  general  set  up  established  in  the  previous  paragraphs  leads  to  a)  the 
definition  of  a boundary  value  problems  for  compatible  elements,  associated 
with  the  given  Gauss' formula;  b)  the  possibility  of  deducing  classes  of 
variational  formulations  of  this  boundary  value  problem. 


By  definition,  a compatible  element  £ r C is  uniquely  characteri- 
zed by  a constraint  element  C of  the  type  : 


c.  L ifps 

The  element  Z can  thus  be  interpreted  as  the  (generalized)  solution 
of  the  following  boundary  value  problem  ; 


6 


i 


\ 


O //) 


- 6 a - 


41 

whereV  denotes  the  set  generated  by  C.The  notation  n(s)  is  meant  to 
indicate  that  for  problem  (1)  the  values  of  I'on  7)^.  J2.  are  uniquely 

determined  once  s is  given.  ^ 

The  corresponding  domains  are  (closed)  subspaces  of  the  solution  space 
referred  to  as  basic  subspaces  and  denoted  by  : 

4 , F(CJ  ; Xr  - F(Cr)  ■ - F(C*);  4-  F(Ct)  <3*. 

The  isomorphism  F makes  it  possible  to  define  a pseudo-dimensionality 
(psd)  for  subspaces  of  H as  the  number  of  constraints  (i.e.  number  of 
components  of  the  constraint  element)  which  can  be  given,  independently, 
arbitrary  values.  Thus  psd  H = 4,  any  of  the  basic  subspaces  has  psd  = 1 
and  any  fixed  element  "p  6 //  has  psd  = 0^2-7. 

The  solution  space  H has  the  direct  sum  decomposition  ; 

H --  4 co  xF  © 4 © 4 

the  basic  subspaces  being  mutually  orthogonal. 


Given  the  property  of  H,  the  only  thing  which  needs  to  be  proven  is  the 
orthogonality  of  the  basic  subspaces.  This  follows  readily  from  the  de- 
finitions (2.^,8)  of  the  subspaces  and  from  the  Gauss  formulap.^  ) rewrit- 
ten as  to  exhibit  explicitely  the  components  of  the  constraint  element  (4©) : 


C*,  Z’> 


< Ft.  Ft'  > 


Pv,  nj 


7 


The  element  t of  the  vector  space  H will  be  referred  to  as  the  solution 
vector  and  the  values  £ , f7,  rz-  of  the  source,  principal  and  natural 

boundary  constraints  will  be  referred  to  as  the  constraints  "of"  the 
solution  vector  'f 

2.,*-!  Linear  varieties  and  two-field  variational  functional  - 

Let  ■% . ye  be  two  arbitrary  closed  subspaces  of  //'</<*/  j 3 

1 ' /O  O ^ ' 

their  orthogonal  complements  and  [3  ^ the  corresponds 

projection  operators  (e.g.  U • T » (P  //  ) . 


Denote  by  the  linear  varieties  obtained  by  translating  and  c/n 

by  r : 


;lating  s/  and  -/ c 


7 - j <5 u I z r y,  ; £■ c ^ J 

■'t  c 

% -•  / ^ 6 h!  TP  -■  ^ ^ e -Jfi  i 


Since  Z C the  projection  on  of  any  element  ^ 6 

satisfies  the  same  constraints  as  the  corresponding  projection  of  the 
solution  vector  Z . Similarly  for  the  elements  ^ 6-  17a  . The  elements 
of  the  varieties  can  thus  be  loosely  characterized  as  satisfying  a cer- 
tain set  (in  particular  null  if  ty-H)  of  constraints  of  the  solution 4.2  1,1^] 


Define-' tfte-  two-f 


- 8 - 


Define  the  two-field  variational  functional  : 


> 


cp> 


£ V«  j 6 ^ 


(2-/2  ) 


As  proven  in  l_  3 ] this  functional  attains  a stationary  value  for  - 27 ^ 


and  'cp>  - ”2^  such  that  : 


^ * ; 


fp  e - c.  ->  -Y*  n ^ 


This  functional  will  be  used  to  obtain,  via  the  Gauss  formula,  different 
classes  of  variational  problems  equivalent  to  the  boundary  value  problem^//) 


3.-  VARIATIONAL  FUNCTIONALS 


The  Gauss's  formula  (2.5)  is  used  to  express  the  two-field  variational 
functional  (Z,  /£)  in  terms  of  the  components  of  the  constraint  vector  -G  - 

By  recalling  that  Z - ^ V = / /°  / v /<■  3 

and  by  letting  : 

\ - <=  F<  - = Fc * --  Ft  V,  ^ J 

with  a similar  notation  for  2^  it  is  readily  shown  that  : 


k7^,7p  )-■  i * > " ^ 


Z“  -V  > f 


r i 


P- 


4 


- 9 - 


and,  by  further  developing  ^ 


[(?'%,  y - £ />X  ^ ^ 

Each  one  of  the  above  four  expressions  leads  to  possible  forms  of  one  or  two- 

field  variational  functionals  which  include,  as  it  will  be  seen  shortly,  the 

classical  primal  and  dual  one  field  functionals.  In  each  case  there  is  a minimum 

number  of  requirements  to  be  imposed  on  the  varieties  If.  and  IX  in  order  to  el i - 

* P 

minate,  as  obviously  needed,  the  terms  containing  unknown  quantities  related  to 
the  solution  element./  i.e.  quantities  other  than  the  components  of  the  given 

->7  - 

constraint  elementCJ.  Such  terms,  underlined  in  equations  [7>l)lr'car\  be  elimina- 
ted only  by  acting  on  the  component  of  the  constraint  elements  and  Ca 
defining  the  linear  varieties. 

Since  ^ and  ;^are  subspaces  translated  by  the  components  of  £ and/or  to 
be  acted  upon  must  be  set  equal  to  the  corresponding  components  of  C.  . This 
automatically*  identifies  also  the  subspaces  ^ and,  upon  equation  (fo,  t£>). 


the  variety  ]f  characterizing  the  eventual  arbitrariness  in  2:^ . 
rements  can  obviously  be  added  to  further  restrict  the  varieties 


Other  requi- 

LTX  and  tfy 


In  equation  ( 3 / ) the  minimum  requirement  criterion  is  satisfied  if  one  takes 
for  then  the  (^contribution  becomes  constant  and  thus  irrelevant. 

Hence  = ff  • 1 f = //  and  (with  obvious  notations)  : 
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Vp  -f 

Xp  : F j L p,  np/  t/i  J ] ' ^ L •'/  f'  v ' ^ 1 " 

= f j ^ t j f*  ^^4^4 


Thus  : ^ - f^; 

the  functional  reads 


no  indeterminacy  in  c exists  and 

P 


\ k = i <z>*-<Tp.  t>  - 4 ( p^»  .) 

Ice  H ■ Tp  -•  6 rTp  ; f)-  l"p  /,  - /'•  ( /«/* 

V * 


This  is  the  most  general  two-field  hybrid  functional  that  can  be  obtained 
since  is  an  arbitrary  (i.e.  unconstrained)  element  of  H and  £.r  is 
subject  only  to  the  principal  constraints  on  o & . ' 


A similar  procedure  applied  to  equation  (3-2/  ) yields  : 

V = t ; H[i  --  » ; w « '4  = // 

Hp  ; /p(£.  ft  ^ '-'f.  - ~f  Xp(t,  X-t 


and  the  variational  functional  reads  : 


No  indeterminacy  in  2*  exists.  One  of  the  elements  is  still  completely 
unconstrained,  the  otner  satisfies  the  source  and  the  natural  boundary 
constraints  only. 
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When  considering  the  other  two  general  expressions  for  K ( Z ) the  term 

Ci  m2-  f P 

.Tl  'T>may  be  made  to  disappear  since  otherwise  essentially  falls  back  into  the 

P/  V ' r'  A - 

cases  already  dealt  with. 


Accordingly,  the  minimum  criterion  applied  to  equation  (33)  leads  to  : 


X. 


'Tp.o  P.f’ 

, tfp.  * Z r 

and  the  variational  functional  reads  : 

1 V I'vr*  JZ  ~ f}  t1  ~ f ~ 


? ) 


m ) 


This  is  a quite  general  two-fields  hybrid  functional  since,  once  again,  no 
constraints  are  imposes  on  2^  and  only  the  principal  boundary  constraints 
are  imposed  on  l>p  . Compared  with  the  functional  (3  S ),  (Tj^is  disappeared  but 
the  source  constraints  of  2^  appears  explicitely.  If  the  restriction 
imposed  on  , the  term  ( ^ ) disappears  and  : 


l\  . ^ » tr&x&x 


The  indeterminacy  in  is  however  irrelevant  since  as  said,  K would  not  con- 
tain any  volume  integral  of  . The  element  Jp  can  be  made  to  disappear 
all  together  by  further  imposing  the  natural  boundary  constraint  to  the  ele- 
ments Z,  . The  single  field  functional  thus  obtained  is  nothino  but  the  classi- 
cal "dual"  variational  functional  * 


( J 


- i 


N 


- 71  ' 


■Y 


0 vt ' Jl 

) 


i 


(/  v i w ; 


< 
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From  the  last  expression  ( ) one  finally  gets  : 

^ ~ ~ 71 
l ^ - X.j  *,  l p,  - ~ * * p & 

l Zp  f - Z t 


and  the  variational  functional  reads  : 

K-  i <6  ^ ^ S)  - e/>  [p  h,  N'  ^ t'y,  l 


1 


0*4-  ™> 


nj 
C $■  S) 


No  constraint  is  imposed  on  £^and  the  indeterminacy  in  c^is  once  again 
irrelevant.  As  before,  7p  can  be  made  to  disappear  by  imposing  the  prin- 
cipal boundary  constraints  on  lZ  . The  single  field  functional  thus 

P yJ  C . > 

obtained  is  nothing  te£h  the  classical  "primal"  variational  functional  Kp 
for  the  boundary  value  problem  {2,11)  : 


I kf  -.i  <6  v1-  7 ; * Z ( 


v ) - kp/Vu  j 

( 3 ,/£■; 
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9-  GENERALIZATIONS 

The  classes  of  variational  functionals  can  be  ba>  widened  by  further 
decomposing  the  subspaces  Xp  and  • The  argument  will  be  developed 

in  detail  only  for  Xp.  The  constraint  element  corresponding  to  J^is  : 

c P --  L o,  p,  o,  o J 

Let  J"  be  the  index  set  [l,2, ,mjjand  any  partition  of  J and 

put  p = pa  + pb  with  l [>y  j ' r^nd  ^ -- l ty] ; / t V 

The  restriction  of  < y > to  Xp  will  be  written  as  : 

£ ry  K- 7 'J * 5 - *>  ■ ( ?j,  Kf 7 'J 

j -T<.  </  / ' ^ ^ (4./J 

* ck  (K,  ? ( K , Kb  7 ') 


' 'a, 

Each  contribution  on  the  right  hand  side  can  be  interpreted  as  the  restric- 
tion of  the  inner  product  to  two  orthogonal  complementary  subspaces  X« 

r Sc 

h 


and  defined  by  : 


-V.  ■ / r - Xr  I 5'  <•'  h ,■ . a -•  ty  > J e T,  ■ ly  v±  - , h j 

V " 


(Jt-i) 


^s'P/W,  ' fr//6  1 'X 


V V ■ />  v 
^'.•2 


Indeed,  the  corresponding  constraint  elements  read  : 

C K ^ c?/  ^ '’c"  y 6'  ^ 

<->  - Ph 

' l TV 


1 


p 


- 14 


--  Ct 


C 


v*  c £■ 


t> fc  't 


the  decomposition 


so  that  ' ""fit  T Lfb 

being  unique  when  the  usual  assumption  upon  the  existence  and  uniqueness 

of  the  solution  of  the  boundary  value  problem  corresponding  to  X ^ 

i h 


vU  . 


Hence 


and  the  orthogonality  follows  from  equations  (4X  ) and  (V.2)  • Notice 
that  for  any  Z 6 X p it  is  h/.'  25  tr=  C'  (</  6 J^Jon  the  entire  bounda- 

(P 


ry  'y  Jc  since  for  any  Xp  it  is  hi  Z - Y1  = O . Either  subspa- 
ce Xpu  or  Xp  can  be  further  decomposed  with  the  same  procedure  but 
such  a development  will  be  omitted. 


Consider  now  the  functional  K ( r, 
on  the  variety^  Ip,  by  lettijng  : 


) and  weaken  the  constraints 


V 


LX 


5 J ’ '!i 


- 


where  c0is  an  arbitrary  but  fixed  element  of  H. 

Hence  -^pis  translated  by  an  element  ct  whose  projection  on  ^ is  no 
longer  equal  to  that  of  the  solution  element  but  quite  arbitrary.  It  is 
readily  shown  that  the  stationary,  properties  mentioned  in  paragraph 

still  hold  with  V^,  the  translate  of  r\  ^ by  Zt  instead  of  Z. 

As  a consequence  of  this  new  definition  of  LpOne  has  "more  freedom"  in 
the  process  of  elimination  of  the  underlined  terms  from  equations  (i  / )~  (3xJ. 
Indeed  one  can  set  the  components  of  the  constraint  element  ^correspon- 
ding to  the  projection  fp  equal  to  some  fixed  but  otherwise  arbitra- 
ry value.  As  this  can  be  done  only  when  J-  , i.e.r1/^/  ft  , one  rea- 
dily sees  that  the  suggested  procedure  amounts  to  "shift"  a constraint 
from  £pto  . 


L 
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Thus,  for  instance,  one  can  eliminate  the  underlined  term  from  equation 


with  arbitrary  but  fixed.  The  variational  formulation  is  equivalent  to 
the  original  boundary  value  problem  if  one  restricts  to  the  subset 
satisfying  the  principal  boundary  conditions  h . 


Hence  the  new  variational  formulation  reads  : 

K"  - { <•  Op,  T ) ' ( v'irp/  * ) 


Vp  , - r 

o r top  -Ji  V 

•7 


f;  >) 


T 


c-  X 


Comparison  with  equations  (3.5*)  shows  that  the  principal  boundary  constraint 
has  been  split  into(^t  and  the  part  j?  has  been  "shifted"  from  the  ele- 
ments to  the  elements  2”  . I® 

r °< 


Similar  developments  can  be  applied  to  the  other  three  exoressions  (3  £), 
(33),  ( 3 y)  and/or  for  the  natural  boundary  constraints  to  obtain  correspon- 
ding new  sets  of  variational  functionals. 


CONCLUDING  REMARKS 

The  inverse  approach  developed  here  makes  it  possible  to  unitarely  derive  hybrid 
and/or  primal  and  dual  variational  formulations  of  classes  of  boundary  value 
Droblem^by  starting  from  general  explicit  expressions  for  Gauss'formulae.  The 
"degrees  of  freedom"  are  essentially  constituted  by  the  different  "physical 


w 
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nature"  of  the  elements  £ (e.q.  symmetric  second  order  tensors,  general  second 
order  tensors,  hiqher  order  tensors,  and  so  on),  by  the  definition  of  the  inner 
product  (which,  amonq  other  thinqs,  characterizes  the  set  of  principal  and  na- 
tural boundary  differential  operators  and  hence  the  types  of  boundary  conditions 
that  can  be  analysed)  and  by  the  selection  of  the  "constitutive"  part<$lTof  2C 
^which,  for  second  and  hiqher  order  tensors,  can  be  done  in  more  than  one  way j . 
Chanqinq  any  one  of  the  above  "deqrees  of  freedom"  leads  to  different  classes 
of  boundary  value  problems. 
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THE  CLOSED  SPLINE  FUNCTIONS 

by 

Luigi  G.  Napolitano  - V.  Losito 


1.-  INTRODUCTION 


The  interpolation  problem  solved  by  the  classical  spline  functions  is  well 
known  [l  ] . 

Given  (n)  points  (t^) 

a < tj.  < t* < t/n  * k 

in  the  closed  interval  j^a,  b]  and  (n)  real  numbers  (r^)  the  interpolating 
spline-function  of  order  q£  n corresponding  to  the  (n)  points  (t.)  is 
the  only  function  which  solves  the  following  minimization  problem  : 

b (l)  - 

/rvww,  ■ ■ ' 


-f 


[4  («]  M 

«. 


1 £ l £r 


m 


(1.1) 


where  h’M  is  the  Hilbert  space  of  real  functions,  defined  on  0 * bJ 
having  a square-integrable  q-th  derivative.  The  space  Hq  quantifies  the  notion 
of  "degree  of  smoothness"  of  the  interpolation  curve.  For  q = n the  interpola- 
ting spline  reduces  to  the  unique  polynomial  of  degree  (n)  passing  through  the 
points 

The  interpolating  spline  belongs  to  the  subspace  of  real  functions  s(t) 


a)  Airfoil 
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b)  Spline-scheme 

defined  on  fa,  bj  and  such  that  f'3: 

a)  s is  a polynomial  of  degree  (2q  - 1)  in  each  of  the  open  intervals 
( 1 £ i 6 n - 1 ) 

b)  s is  a polynomial  of  degree  (q  - 1)  infl^f  andjt  0[y  i kl 

c ) s^q  ^ is  continuous  on  , bj 

In  many  aerospace  problems  it  is  necessary  to  find  interpolating  functions 
for  "airfoils",  i.e.  for  closed  curves.  This  interpolation  problem  is  often 
solved  by  using  the  classical  spline  function  with  fcj = (see  figure). 

Such  a procedure  is  unsatisfactory  for  two  reasons.  Firstly,  the  discontinui- 
ties introduced  at  the  airfoil  point  (usually,  the  trailing  edge)  correspon- 
ding to  the  values  and  t^are  not  those  proper  of  the  airfoil  (if  any). 
Indeed  the  classical  spline  of  order  q yields  different  and  uniquely  defi 
ned  values  for  the  first  (q  - 1)  derivatives  at  the  points  t = t-j  and  t = tp. 

Secondly,  upon  the  property  b)  of  the  spline  functions,  all  derivatives  from 
order  q to  2 (q  - 1)  of  the  interpolating  function  will  vanish  at  the  point.  P. 
This  expecially  for  low  values  of  q,  much  to  restrictive  in  practical  applica 
tions . 

In  principle,  this  approach  could  be  improved  upon  by  imposing  a number  of  addi 
tional  requirements  at  the  points  t = t1  and  t = t to  better  approximate  the 
shape  of  the  airfoil  at  the  points  P.  Not  much  could  however  be  said  about  the 
properties  (such  as  existence,  uniqueness,  extremality,  and  so  on)  of  the  re- 
sulting classes  of  interpolating  functions. 
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A completely  different  approach  is  advocated  here.  Rather  than  first  assuming 
that  the  points  t ' belong  to  an  interval^ t>J  and  then,  somehow  accounting  f< 
the  fact  that  both  and  correspond  to  a same  point  of  the  curve  to  ' 
approximated,  the  points  ( ) are  right  from  the  beginning  suonosed  to  be 

to  an  (arbitrary)  closed  contour  C. 

By  thus  doing  each  t-  corresponds  to  one  and  only  one  noint  of  the  curve  to  r» 
approximated  and  the  interpolating  problem  for  closed  curves  is  conseouentlv 
cost  into  its  more  natural  formulation. 

If,  the  present  paper  it  will  be  shov/n  that  the  interpolation  problem  thus  for-” 
lated  is  amenable  to  rigorous  analysis  and  leads  to  a new  class  of  spline-fur 
tions  defined  on  a closed  contour  C. 

In  wanting  of  a better  terminology,  they  will  be  referred  to  as  "closer  , 
functions"  since  they  yield  interpolating  functions  for  closed  curves 

Hilbert  space  theory  of  spline-functions  will  be  used.  Basic  results  needed 
recalled  in  paragraph  (2).  Existence,  uniqueness  and  characterization  of  the 
closed  spline  functions  are  established  in  paragraph  (3). 

These  results  are  restated  in  conventional  (i.e.  non  abstract  snace)  notation 
in  paragraph  (4)  where  additional  relevant  properties  (following  from  the  ha;  • 
abstract  theory)  are  also  given.  A final  paragraph  offers  few  comments  on  the 
other  classes  of  closed-spline  functions  which  can  be  constructed  by  the  same 
approach . 
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2.  - THE  ABSTRACT  SPLINE  FUNCTION  THEORY  - 

The  needed  basic  results  of  the  Hilbert  space  formulation  of  spline-function 
theory  are  summarized  here  for  ready  and  convenient  reference  £ See  [l] for 
greater  detailsj  . 

Let  X,y,  be  three  real  Hilbert  spaces  with  norms 

and  I ’•  X Y j A ; X % two  linear  continuous 

operators,  which,  without  any  loss  of  generality  are  supposed  to  be  o*t  • 

Given  any  fixed  *e 2 , define  the  set,  supposed  non  empty: 

I?  ' t X / c (2.1) 

and  consider  the  following  minimum  problem  : 

1/ Tally  - 

The  element  S'fc  C ^ , if  it  exists,  is  called  the  interpolating  spline 
corresponding  to  ( T/  ). 

The  following  existence  and  uniqueness  theorem  holds  W&onom  ♦ ( Exi'a'tenee) . 
Theorem  1 . 

"The  solution  of  problem  (2.2)  exists  for  any  < t iff  N(T)  + N (A)  is 
closed  in  X and  is  unique  iff,  in  addition,  N(T)  AN(A)  = ^0^ 


1 


Here  N(B)  denotes  the  null  space  of  the  operator  B. 
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The  definition  of  the  space  S of  spline-functions  follows  from  the  characteriza- 
tion stated  in  : 

Theorem  2 : 

"Under  the  hypothesis  of  Theorem  1 , (j  £ is  the  interpolating  spline  cor- 
responding to  (T,M  ) iff 

^(5-T*  >y  = O v Xfc  A/W)-I0 

Here  j >y  denotes  the  inner  product  in  V . Similar  notation  will  be  adop- 

ted for  the  inner  products^. 


The  space  S is  the  subspace  of  X defined  by  : 

5-loet  I <T*,TxVo  ;Vxe/VfA)j 

Theorem  1 can  be  formulated  (and  generalized)  in  S as  : 


(2.3) 


Theorem  3 


"For  any  fc  ^ there  exists  a unique  element  6*6  5 such  that  A(«J  = z". 

The  extremal  properties  of  interpolating  spline  are  condensed  in  the  following 
theorem. 


Theorem  4 

If  z is  an  arbitrary  fixed  element  of  Z and  (T  is  the  unique  element  of  S such 
that  A ( 6"  ) = z then  : 
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jT(f-x)|)^  ^ \\rb-*)\\y 

V x t Ijj 

and  any  other  5-^  S having  this  property  belongs  to  the  set 

[?}=  <t+N(t) 

||  t («•- <>)||  = ||t(x-«)|| 

11  y xtit  y 

and  <T  is  the  unique  element  of  1,  having  this  property. 


(2.4) 


(2.5) 


Further  discussion  will  be  limited  to  the  case  of  more  specific  interest  here, 
in  which  the  number  of  interpolating  constraints  is  finite. 

The  space  Z is  consequently  finite  and,  upon  theorem  1,  existence  is  automati- 
cally guaranteed  whereas  uniqueness  requires  that  N(T)  be  also  finite. 

'W 

Suppose  then  that  Z t R ,with  the  usual  inner  product,  and  that  N(T)  is  of 
dimension  q.  Then  : a)  the  operator  A can  be  expressed  as  : 


Ax  = K±,*>x  >•••••/<  >*^6  f? 


(2.6) 


where  then's  are  n independent,  linear  continuous  functional  on  X;  b)  the  set 
can  be  written  as  : 


x t y /<  <lyx>x  = 


= *1  ) J * i-s  v*  fa, R 


c)  the  following  characterization  theorem  holds  : 


Theorem  5 : 


" is  the  interpolating  spline-function  corresponding  to  (Tyk^l;  ) 

iff  there  exist  (n)  coefficients  X • such  that  : 


(2.7) 
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Here  T is  the  adjoint  of  T and£  U denotes  the  orthogonal  subspace. 

The  different  classes  of  spline  functions  that  can  be  obtained  from  the  above 
abstract  formulation  depends  on  the  choices  of  X,Y  , T and  K;  . 

A particular  instance  is  developed  in  the  next  paragraph. 

3 - The  closed  spline-functions  - Existence,  uniqueness  and  characterization 


Let  C be  a closed  contour,  sufficiently  smooth  and  regular  for  the  validity  of 
all  properties  that  will  be  used. 

°l 

Take  X = H (C)  and  Y = H°(C)  with  their  standard  inner  product. 


Denote  by  "S'  the  curvilinear  coordinate  along  C measured  from  an  arbitrary 
initial  point  P0  and  normalized  with  respect  to  the  length-^  of  C. 


Thus  0***1  and  't  (P0  ) r 0*  j )- 

As  operator  T:  take  ol  where  (d)  stands  for  the 

derivative  with  respect  to  “Zi  . Clearly  T is  linear  and  continuous. 


Finally,  for  Z c K the  (n)  functionals  : H 


be  defined  by  : 


(3.1) 


0 = 


i 


with 
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The  null  spaces  of  the  operators  A and  I are  then  given  by  : 

Hl(c)  j 

N(t)*  J X t H\c)J  ojs|*6  ti\c))x.--  C^Jj  (3.2) 

Hence  their  dimensions  are  equal  to  (n)  and  to  one,  respectively,  and 
N(r)n  P((fi)  ~ jot  provided  n > 1 . 


The  minimum  problem  (2.2)  reads  : 


- (*ti  f.  ...  ,Kf^)  (:  R 

and,  according  to  theorem  1,  it  has  a unique  solution,  for  any  r, 
n >1. 


(3.3) 


as  long  as 


The  existence  and  uniqueness  of  s(v)  having  been  established,  the  next  basic 
task  is  to  characterize  it. 


Upon  the  definition  of  adjoint,  theorem  5 and  equations  (2.7),  (3.1)  one 
has,  subsequently  : 


Cj  Cj  / 1 1 

< d <r,  d * >H„  = < (d  J d r,  * = 

£ nv 

r Z Z,  K;  jx  = 2 ^ x fat ) v x t fi  (c) 


(3.4) 
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Always  upon  theorem  5,  l ^ L 

ir'/ 

so  that  the  coefficients  /-must  be  such  that  : 

c 2.  /t  = 0 J 

l- 1 n 

Now,  as  seen,  oU  a/CtJ-  d so  that  always  upon  equation  (3.2)^  it 
suffices  to  impose  that  ^ be  ff*i  orthogonal  to  unity. 

On  account  of  eq.  (3.1)  this  leads  to 


• 'Is 

2^ 


(3.5) 


as  the  only  homogeneous  condition  to  be  satisfied  by  the  X't  .^This  is  to  be 
contrasted  with  the  situation  of  classical  spline  - functions  for  which  the 
coefficients  must  satisfy  the  (q)  homogeneous  conditions  ^ /(•  C 

[os.  K<  ej-i)] 

In  order  to  characterize  CT^  or,  equivalently,  to  obtain  the  additional  con- 
ditions which  defines  the  4 uniquely  one  must  express  the  rightmost  term 
in  equation  (3.4)  as  a scalar  product  in  H°(C)  valid  for  any  Xk  H\c). 

A Mac-Laurin  development  ofx(^J  with  the  rest  expressed  in  intergral  forms 


where 


+ i (?c-yl  x%)jt 

j <H  J h-w 

(4  i -Izo 

(*  ) * T 

* ) <(  i-o 


0 
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Thus,  on  account  of  eq.  (3.5)  : 

* 


2 k *(*.•)- 


,^-f  1 X^o)  ^ ^ (3.6) 

> I “*  J‘.  fr.  l \ / 


Cf-U' 


‘V 

Upon  the  (absolute)  continuity  of  the  derivatives  * (?)  for  1 £ k£  q-1 
one  can  let  : 

1 2. [t  *>■  I'*J *%) - x^t) a* 

‘ )■*.  i • (3.7) 

2.  h- 

»*'  . 
where  the  (q-1)  coefficients  V_  are  uniquely  determined  in  terms  of  A^, 

and  ^ . It  is  convenient  to  express  this  functional  dependence  in  terms  of 

the  discontinuities  of  'iTf'jrj  and  of  its  first  (q  - 1)  derivatives  at  the 

"initial"  point  f0  , given  by  : 

'/K}r)-rw(J*)=  \ h cfk<a-i  (3.8) 

j**«  trio! 

Repeated  integrations  by  parts  lead  to  : 

Co)  y x.  e HVC) 

J K^o 

so  that,  by  comparison  with  the  first  of  equations  (3.7)  : 

/ . \K  c C *)  ^ 

('<)  ot  = / AlZL—  o^Ksq-i  <3-9) 

i-T  (VK-l)l 

r ft-0 

StT  - C> 
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The  last  of  these  equations  is  identically  satisfied,  upon  equation  (3.8) 
and  the  others  yield  the  required  (q  - 1)  equations  needed  to  determine^  ir, 
terms  of  and  ^ . Equations  (3.7)  and  (3.9)  affords  the  required  transfor- 
mation for  the  first  term  on  the  right  hand  side  of  eq.  (3.6).  The  transforma 
tion  of  the  second  term  hinges  on  the  identity  : 

fo-'fe)’  *=  (-o1  fc"*? 1 

Let,  on  account  of  eq.  (3.5)  : 


■ (i-O!  i‘°  jj! 


(3.10) 


with  : 


d ; = (-0  2.  1 

w h-1-0! 


(3.11) 


Since  : 


Co ) V*  i U (c),v 


(3.12) 


The  term  ino(^does  not  contribute  to  the  integral  appearing  in  equation  (3.6) 
and  one  may  thus  let  : 

1=1  (q~0!  fro!  r 


K*)<-0 


r r- 
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Substitution  of  equations  (3.13),  (3.7)  into  equation  (3.6)  leads  finally 
to  : 


2 A *(?<) / v x b ^ ^ 

y 4’)1  2 £(*£$'+!>('&)+**) 


(3.14) 


Comparison  between  eqs.  (3.14)  and  (3.4)  shows  , on  account  of  eq.  (3.12), 
that  . \ /»  Q H ' 


tnat  q . / \7 

A)=  <f+(W-o  l , j&s)  (?_i5) 

a } 1 

where  p is  an  arbitrary  constant.  « 

fi) 

The  expression  (3.  15)  for  « contains  (n  + q)  coefficient  (n  coefficients  A' 
and  q coefficients  ).  The  following  lemma  shows  that,  on  account 

also  of  equation  (3.5),  only  (n)  of  them  are  arbitrary. 

Lemma 

"The  function  (T  ('fc)  and  its  first  (q  - 2)  derivatives  are  continuous  on  C" 
Proof 

Given  the  properties  of  polynomials  and  of  the  functions  ("£  - ^ ) the 
statement  of  the  lemma  needs  to  be  proven  only  at  the  point  Po. 

In  the  interval  Dm  o,  \ J eq.  (3.15)  reduces  to  : 

* b A s!!'  + IL  + ifr)+r(t)  > r e[o,xJ  <316> 
(l-O!  1 
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‘n  the  interval  yif  fig  - for  any  i so  that  : 

(*l)  °\ 

g*  (%)=  ^ + (- 0 ^ tf&)  ; t e (3.i7) 


S'r.e  for  ('£- 'fcjj | _ (T-'fci')  upon  eqs.  (3.10)  and  the  second  of  eqs. 
(3.13)  : 

< f 2 +*b>S  d % r (-  I )\  - p(  V) 

M1  j=»  j.'  ? 

From  eqs.  (3.16)  and  (3.17)  it  follows  that  : 

. (ykl  0*“'  , fK)  it)  . , M <i 

acr  (o*j-<r  (c)z  £V _j/  6»)«  <5ir  +(-<)<,„< 

k t[o,vlI 

since,  upon  the  second  of  eqs.  (3.13)  k(o)=.  0 and 

Hence,  upon  equations  (3.9)  and  (3.11)  : 
and  the  proposition  is  proved. 


rhe(q  - 1)  conditions 


r (V*)  - 

6<T  *•  O 
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yield  the  (q  - 1)  equations  which,  together  with  equation  (3.5),  make  it 
possible  to  eliminate  q coefficients 
which  depends. 

J1) 

The  (q  - 1)  - th  derivative  of  O is  constant  in  each  open  interval/'?'  % j 
1 £-  i < n).  Its  discontinuity  at  the  points  of  coordinate  ?■  is  equal  to  / . 

" l 


( 


) out  of  the  (n  + q)  upon 


The  closed  interpolating  spline  <5"  ( ) is  obtained  by  integrating  equa- 

tion (3.15)  q times.  The  corresponding  (q)  arbitrary  constants  are  determined 
by  imposing  the  continuity  of  (S'  and  of  its  first  (q  - 1)  derivatives  at  . 

The  function  thus  obtained  will  contain  (n)  arbitrary  constants  which  are 
determined  (uniquely,  upon  theorem  5)  from  the  (n)  conditions  Gife)  = ’tj, . 

Existence,  uniqueness  and  main  characterization  of  the  closed  spline-functions 
have  thus  been  proved. 


In  the  next  paragraph  the  additional  properties  of  these  classes  of  functions, 
as  derived  from  the  general  theory  outlined  in  paragraph  (2),  will  be  pre- 
sented. To  facilitate  a more  wide  spread  comprehension  abstract  space  nota- 
tion will  be  abandoned.  To  give  a complete  panorama,  the  initial  interpolation 
problem  and  its  solution  will  also  be  reformulated  in  standard  notation. 


4 - The  Closed  spline  functions  - Definition,  characterization  and  properties 


Given  a sufficiently  smooth  closed  contour  C and  (n)  arbitrary  points  P. 
on  it,  let  ( i £ l £ m. J be  their  curvilinear  coordinates  (see 

figure)  measured  from  an  arbitrary 
initial  point  P and  normalized  with 
respect  to  the  length  of  C. 
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Define  the  (n)  open  intervals  ,r 

and  prescribe  (n)  values  *£•  with 

fo,**, 

The  closed  interpolating  spline  function  S'"  of  degree  (q)  corresponding 

t0  1 \ £ i £ m ) is  defined  as  the  unique  element  of^(C)  such 

that  : 


•^J^s  - J[f ?!]  ^ 


It  = i hUc)/  li  > ( 

In  order  that  <T  t H (cj  be  the  closed  interpolating  spline-function 
corresponding  to  ( ) it  is  necessary  and  sufficient  that  : 


a) <o  be  a polynomial  of  degree  (2q  - 1)  in  each  open  interval  7 fcV  X 

(Vo;V.-iJ  y L 

b)  <T  be  continuous,  on  C,  together  with  its  first  2(q  - 1)  derivatives;  thus, 
in  particular  : 

Xh:)  =X(t.-)  ; VKe[o,Z)-iJ  ; V Zc 

c)  ^ be  such  that 


<*•  ft;  ) = \ 


V c e [i,  -nj 
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The  set  of  functions 
a subspace  S (the 
of  dimension  n. 


^ 6 fi  satisfying  the  first  two  conditions  constitute 
space  of  closed  spline  functions  corresponding  to  Z"' 


) 


For  any  element  'l*-  S (cfr.  theorem  2)  : 

A base  of  S is  given  by  (n)  A-  4 such  that  (t  ) =.  £ and  it  is  : 

^ ^ ^ 
i=i 

The  following  extremal  properties  hold  (cfr.  theorem  4)  : 


Given  (n)  arbitrary  but  fixed  , if  <T  is  the  unique  element  of  S such 

that (T ( ) = Y-  then  : 

* v 


A) 


[c%) 


Z fr-  (<()  Jl) -Z 


t)  Id'z-  (tajm 
J 4 fe  5 


and  any  other  element  Ofe  S having  this  property  differs  from  G by  a 
constant  : 


B)  for  any  /Jfc$  : 


-£j 


and  6"  is  the  unique  element  of  having  this  property. 

An  element  4 6 HYO  belongs  to  the  space  S of  the  closed  spline-functions 
corresponding  to  2"^  iff  it  is  representable  as  : 

j-o  }! 

where  the  (n)  coefficients  and  the  (2q)  coefficient  ^ satisfy  the  follo- 
wing 2q  equations  : 


2A- 


Zfl-1  _ -n.  . 2<j  - J 7 

5 A'  ^ I- 0 

^x(rK)j  pit 

05  K < 


where  : 


= A%')  - 

For  any  function  •feHVc)  it  is  : 

'^k) /%)<*?  = ZAffe) 


■C  / 

The  (n)  coefficients  ( Aj ) represent  the  values  of  the  discontinuity  of 
(2  q - 1 ) -th  derivative  of  s at  the  points  of  coordinate 

/ 
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When  the  following  (n)  equations  : 

*(%)**{  (imto) 

are  added  to  eqs.  (4-1),  the  resulting  system  of  (n  + 2q)  equations,  linear 
in  /^and  ^ , admits  an  unique  solution. 

It  may  be  instructive  to  compare  the  closed-spline  functions  considered  herein 
with  the  classical  spline  defined  over  the  interval  fa,  bJ  = [o,lJ  and  corre- 
sponding to  the  partition  : 

0 - A -=ts  ■ 

In  the  classical  spline-functions  the  coefficients  & are  identically  zero  for 
^ 2(j  i anc*  t^ie  following  (g)  relations  nold  f 1 J : 


0 S 


* K < cj~L 


Upon  letting  : 


(SI 


) and  ( 


one  deduces  from  equations  ( ) , ( 

1-1 


) that 


l=K+i 


oiks  2o-Z 
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so  that,  upon  noticing  that  : 

‘ j - / \o)  * < 


& 


one  finds  that 


7 00 


Ctc), 


Cj  <.  )<  i-  Z(j'  Z 
0 5 k £ / 


6 i 0 

with  the  ^ 'uniquely  determined  for  0*  ^^^-I.Thus  one  recovers  that, 
as  said  in  the  introduction  if  one  simply  takes  and  imposes 

r Km  = -d  the  first  (q  - 1)  derivatives  of  s ( % ) are  discon- 
tinuous at  and  the  other  derivatives,  up  to  the  order  2q  - 2,  vanish 

identically. 


- CONCLUDING  REMARKS 

As  with  the  classical  spline-functions,  other  classes  of  closed-spline  functions 
can  be  introduced  by  making  different  assumptions  as  to  the  operator  A.  Thus, 
for  instance,  one  could  construct  Hermitian  closed-spl ines  £by  imposing, 
as  interpolating  constraints  at  the  points  , also  the  values  of  the  first 
derivative  of  the  function  If  (.c)]  mixed  closed-splines  constraints 
equal  to  linear  combinations  of  values  for  and  ^ , Fourier  closed-spl ined , 
and  so  on.  Furthermore,  one  may  construct  similar  classes  of  interpolating 
and  approximating  "closed-splines  (cfr.  Til). 

These  further  developments  will  be  presented  in  future  reports. 


